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Abstract 

Many fruit trees undergo cycles of heavy fruit load (ON-Crop) in one year, followed by low fruit load (OFF-Crop) the 
following year, a phenomenon known as alternate bearing (AB). The mechanism by which fruit load affects flowering 
induction during the following year (return bloom) is still unclear. Although not proven, it is commonly accepted that 
the fruit or an organ which senses fruit presence generates an inhibitory signal that moves into the bud and inhibits 
apical meristem transition. Indeed, fruit removal from ON-Crop trees (de-fruiting) induces return bloom. Identification 
of regulatory or metabolic processes modified in the bud in association with altered fruit load might shed light on the 
nature of the AB signalling process. The bud transcriptome of de-fruited citrus trees was compared with those of ON- 
and OFF-Crop trees. Fruit removal resulted in relatively rapid changes in global gene expression, including induction 
of photosynthetic genes and proteins. Altered regulatory mechanisms included abscisic acid (ABA) metabolism and 
auxin polar transport. Genes of ABA biosynthesis were induced; however, hormone analyses showed that the ABA 
level was reduced in OFF-Crop buds and in buds shortly following fruit removal. Additionally, genes associated with 
Ca 2+ -dependent auxin polar transport were remarkably induced in buds of OFF-Crop and de-fruited trees. Hormone 
analyses showed that auxin levels were reduced in these buds as compared with ON-Crop buds. In view of the auxin 
transport autoinhibition theory, the possibility that auxin distribution plays a role in determining bud fate is discussed. 
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Introduction 

Fruit trees exhibit two major multiannual reproductive strat- 
egies (Goldschmidt, 2013). In the first, the amount of fruit 
produced allows a sufficient amount of vegetative growth to 
support production of an ample number of flowers during the 
following year (return bloom). Such trees, including fig and 
some orange and grapefruit cultivars, are defined as regular 
bearers. They are characterized by a relatively stable multian- 
nual yield, and usually possess efficient mechanism(s) to con- 
trol excess fruit production. A second strategy is also used by 
trees that bear a heavy fruit load (ON-Crop) in one year, which 
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inhibits return bloom and vegetative growth the next year 
(Monselise and Goldschmidt, 1982). Thus, the second year 
is characterized by low yield (OFF-Crop) and high vegeta- 
tive growth. Such trees, including olive, pistachio, mandarins, 
and many others, are defined as alternate or biannual bearers 
and they are usually characterized by low self-thinning abil- 
ity (Goldschmidt, 2013). Alternate bearer cultivars present a 
serious economic problem to fruit growers. Therefore, chemi- 
cal or manual fruit thinning are common practices in their 
cultivation (Dennis, 2000). In citrus culture, low temperatures 
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during the autumn and winter are a major factor in inducing 
flowering (Valiente and Albrigo, 2004; Knauer et al., 2011). 
Optimal flowering density is achieved only upon accretion 
of sufficient cool hours. It is assumed that a heavy fruit load 
prevents recognition of the low-temperature flowering induc- 
tive signal and/or blocks later stages of inflorescence, such as 
bud break (Albrigo and Galan-Saiico, 2004; Verreynne and 
Lovatt, 2009). As expected, fruit load affects the expression 
of flowering control genes, FT, LFY, API, TFL, and miR156- 
regulated SQUAMOSA PROMOTER BINDING (SPL5) in 
leaves and buds of citrus (Munoz-Fambuena et al, 2011, 
2012ft; Shalom et al., 2012) as well as in mango and apple 
(Kotoda et al, 2010; Nakagawa et al., 2012). 

The mechanism by which heavy crop load affects return 
bloom is not fully understood. The developing fruit provides 
a strong sink for photoassimilates. It was therefore thought 
that depletion of photoassimilates, especially carbohydrates 
from the bud, prevents flowering induction, a hypothesis 
known as the nutritional theory (Goldschmidt et al, 1985; 
Goldschmidt, 1999). Sucrose was shown to play a regulatory 
role in Arabidopsis flowering control (Eriksson et al, 2006), 
but whether sugars indeed play a regulatory role in flower- 
ing induction under various fruit loads in fruit trees has been 
a controversial issue for many years (Hilgeman et al, 1967; 
Jones et al, 1970, 1974; Goldschmidt and Golomb, 1982; Li 
et al. , 2003a, b). Recent work has shown that trehalose metab- 
olism and its product trehalose-6-phosphate were involved 
in flowering control in Arabidopsis (van Dijken et al., 2004; 
Wahl et al, 2013). It was also shown that two genes encoding 
enzymes associated with trehalose metabolism were induced 
in OFF-Crop buds (Shalom et al, 2012). In addition to the 
nutritional control of alternate bearing (AB), it might well 
be that the fruit itself, or an organ which senses fruit pres- 
ence, generates an inhibitory signal (AB signal) which moves 
into the bud and prevents flowering induction (Bower et al. , 
1990; Talon et al, 1997). Fruit thinning or complete removal 
(de-fruiting) from ON-Crop trees induces return bloom 
(Monselise and Goldschmidt, 1981), thus providing support 
for this notion. Gibberellin (GA) is known to inhibit flow- 
ering in many perennials (Goldschmidt and Samach, 2004; 
Bangerth, 2009). However, while exogenous application of 
GA prevents flowering (Goldschmidt et al, 1997; Munoz- 
Fambuena et al, 2012a; Goldberg-Moeller et al, 2013), the 
question of whether GA acts endogenously to inhibit flower- 
ing is still open. The involvement of abscisic acid (ABA) in 
the regulation of return bloom is even less clear (Jones et al. , 
1976; Goldschmidt, 1984; Koshita et al, 1999; Okuda, 2000). 
Polar auxin transport from a dominant sink was also sug- 
gested as a possible mobile signal affecting flowering (Caaejas 
and Bangerth, 1997; Smith and Samach, 2013). 

Fruit load might act at various developmental stages such 
as flowering induction, transition of the shoot apical meris- 
tem, and subsequent stages of flower development and bud 
break (Verreynne and Lovatt, 2009). Regardless of the source 
of the AB signal and its nature, it must be recognized by its 
receptor in the bud which in turn must make the 'decision' 
of whether to proceed to inflorescence or not. In order to 
investigate metabolic and regulatory processes taking place 



in the bud and affected by fruit load, the transcriptome of 
buds from ON- and OFF-Crop trees was recently compared 
during three developmental stages. Changes in metabolic and 
regulatory pathways, including photosynthesis, and in flavo- 
noid and trehalose metabolism were identified (Shalom et al. , 
2012). However, this work was biased due to the use of an 
Affymetrix Citrus Gene-Chip array that contained ~15 500 
genes. In fact, with the exception of trehalose metabolism, 
no other regulatory pathways were identified. In the current 
work, a complementary approach was taken by comparing 
the transcriptome of buds of de-fruited trees with those of 
ON-Crop trees. The genomic analysis was non-biased, as it 
was based on RNA-deep sequencing. It was possible to iden- 
tify an increase in ABA-metabolizing genes, accompanied by 
a decrease in ABA levels and those of its catabolites in buds 
of de-fruited trees. Moreover, a remarkable increase in the 
expression of genes encoding proteins associated with cal- 
cium-dependent auxin polar transport and a reduction in bud 
endogenous auxin levels following de-fruiting were identified. 
The results are discussed in light of the previously suggested 
auxin transport autoinhibition (ATA) theory (Bangerth, 
1989) and its role in AB (Smith and Samach, 2013). 



Materials and methods 

Plant material and sample collection 

Plant material was collected from a commercial orchard of 1 5-year- 
old Murcott mandarin (Citrus reticulate Blanco) trees grafted on 
sour orange (Citrus aurantium L.), located in the central coastal area 
of Israel, during the years 2011 (an ON year) and 2012 (an OFF 
year). Although most of the trees in the orchard bore similar yields 
in a given year, some were aberrant and showed an opposite AB 
trend. These and nearby trees with the opposite yield status were 
selected. Overall, nine triplets of trees were chosen, with each triplet 
(two ON trees and a nearby OFF tree) being considered one bio- 
logical replicate. Fruits were completely removed (de-fruiting) on 22 
August 201 1 from one of the ON trees in each triplet, and this tree 
was labelled DEF. Samples were collected 1 d prior to de-fruiting 
(Time 0, from ON and OFF trees), 1 week following de-fruiting 
(Time 1, from ON and DEF trees), 2 weeks following de-fruiting 
(Time 2, from ON and DEF trees), and 4 weeks following de-fruit- 
ing (Time 4, from ON, OFF, and DEF trees). 

The three most extreme conditions of spring flush flowering were 
compared (see fig. S2 in Shalom et al, 2012): fruit-bearing flush of 
an ON tree, fruitless flush of an OFF tree, and de-fruited flush of 
a DEF tree. Branches of each of these conditions were collected 
from the southeast side of the trees, taken to the laboratory on ice, 
and buds were separated and immediately frozen in liquid nitrogen. 
The percentages of generative, mixed, and vegetative shoots were 
determined for all branches splitting off from one major 50-60 mm 
diameter branch located on the southeast side of the tree, during 
peak blossoming time (usually the first 10 d of April of the consecu- 
tive year). Selection of the branches to be sampled was done prior 
to bud break. 



RNA extraction, quantification, and qPCR analyses 

Total RNA was extracted, treated, and analysed from ~0.2g of fro- 
zen bud tissue, and cDNA was synthesized, as previously described 
(Shalom et al, 2012). Primers for the genes CiFT2, CsLFY, SPL5, 
RbcS, LHCB3, PRK, PSB28, PSAD, SHM, Fd, NCED3, CAX, 
PBPl-likel, PBPl-like2, NPH3, CA-binding EF handl, B-ACTIN, 
and a dual-labelled probe for CiFT2 were designed based on genomic 
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and expressed sequence tag (EST) sequences (Phytozome, http://www. 
phytozome.net/; HarvEST, http://harvest.ucr.edu/) using Primer 3 
software (Supplementary Table SI available at JXB online). Real- 
time PCR was carried out as described (Shalom et al. , 2012). For the 
CiFT2 dual-labelled probe reaction, real-time PCR was carried out 
as described by Shalom et al. (2012). The mRNA levels of trehalose 
biosynthetic genes and flavonoid biosynthetic genes (Supplementary 
Table S2) were determined by nCounter analysis (Nanostring 
Technologies, Seattle, WA, USA) at the VIB MicroArrays Facility 
(Leuven, Belgium), as described by Shalom et al. (2012). 



RNA deep sequencing 

Extracted RNA integrity was determined by Agilent Bioanalyzer 
(Santa Clara, CA, USA) according to the manufacturer's instruc- 
tions. A 2 (j.g aliquot of total RNA from each sample was prepared 
and used for cDNA library constructions using the TruSeq mRNA 
sample preparation kit according to the manufacturer's protocol 
(Illumina Inc., San Diego, CA, USA, REF 15025062). The librar- 
ies (10 pmol) were run on a single read 100 nucleotide run on the 
HiSeq 2000 (Illumina Inc.) on six lanes. Raw fastq files were quality 
checked using FastQC (http://www.bioinformatics.babraham.ac.uk/ 
projects/fastqc), and adaptor sequences were removed using fastq- 
mcf (https://code.google.eom/p/ea-utils/wiki/FastqMcf). They were 
then aligned to the orange (Citrus sinensis) genome database (Xu 
et al., 2013; http://citrus.hzau.edu.cn/orange/) using TopHat2 (Kim 
et al., 2013). FPKM (fragments per kilobase of transcript per mil- 
lion mapped reads) values, which normalize the read count by the 
length of the fragment and the total number of mapped reads, were 
calculated, and differential expression was checked using Cufflinks 
(Trapnell et al., 2010). A hierarchical clustering heatmap and a 2D 
principle component analysis (2D PC A) plot were generated by 
MATLAB (Mathworks, Cambridge, UK), using the log-FPKM val- 
ues of each gene. For gene ontology (GO) analysis, sequences were 
blasted against the UniRef90 database (Suzek et al. , 2007), by which 
a number of GO annotations were derived for each gene. Singular 
enrichment analysis (SEA), which lists enriched GO terms, was per- 
formed [false discovery rate (FDR) <0.05] for the five best GO terms 
of each gene using the AgriGo interface (http://bioinfo.cau.edu.cn/ 
agriGO/index.php). Differentially expressed genes were also func- 
tionally annotated via the automated Mercator pipeline (Lohse et al. , 
2013) (http://www.gabipd.org/biotools/mercator/) and displayed on 
diagrams of metabolic and other processes using MapMan (Usadel 
et a!., 2009; http://www.gabipd.org/projects/MapMan/). 



Hormonal content analysis 

About 50 mg of bud tissue from each tested sample were lyophilized 
and homogenized to a fine powder in liquid nitrogen using a mor- 
tar and pestle. Quantification of ABA, ABA metabolites, and indole 
acetic acid (IAA) was conducted at the National Research Council 
of Canada (Saskatoon, Saskatchewan, Canada) according to pub- 
lished protocols (http://www.nrc-cnrc.gc.ca). 

Protein extraction and western blot analysis 

Proteins were extracted from buds, quantified, and analysed by west- 
ern blot analysis using specific primary antibodies, raised against the 
following proteins: Rubisco complex large (RbcL) and small (RbcS) 
subunits, ferredoxin (Fd), chlorophyll alb protein [light-harvesting 
chlorophyll alb complex II (LHCII)], and Dl protein (PsbA), as 
described by Maayan et al. (2008). 



Statistical analysis 

The statistical analyses used for qPCR results, hormone analy- 
ses, and inflorescences numbers were one-way analysis of variance 
(ANOVA) with Tukey-Kramer multiple comparison tests as imple- 
mented in the software JMP version 10 (SAS Institute). 



Results 

Fruit removal (de-fruiting) induces back flowering 

Normally, fruit load status is similar among most trees 
of the orchard — trees bear either a heavy crop (ON-Crop 
year) or a low crop (OFF-Crop year). A few trees, however, 
show the opposite trend, allowing the collection of samples 
from nearby trees bearing either high fruit load or low fruit 
load. In order to detect changes which might play a role in 
converting ON- to OFF-Crop buds, fruits were completely 
removed from ON-Crop trees in August and buds were col- 
lected 1, 2, and 4 weeks following de-fruiting. The effect of 
the treatment was verified the following spring by counting 
the number of inflorescences and vegetative shoots (Fig. 1). 
Citrus trees bear two major types of inflorescences: genera- 
tive (leafless) and mixed (leafy; flowers and leaves in vari- 
ous ratios). ON-Crop trees had significantly less generative 
inflorescences compared with OFF-Crop trees and DEF 
trees (22% versus 70-75%) and more vegetative shoots (43% 
versus 3-4%). As expected, no significant differences were 
detected in mixed-type inflorescences. Fruit counting during 
harvest time showed that ON-Crop trees yielded 1635 ±98 
fruits per tree while OFF-Crop trees yielded 36 ± 12 fruits 
per tree. 

Fruit removal alters expression of flowering control and 
other genes in the bud 

In order to determine how quickly buds responded to de- 
fruiting, the mRNAs levels of several flowering control genes 
were quantified during the course of the experiment (Fig. 2). 
The Citrus genome contains three FT genes, but only the 
expression of GFT2 correlated well with tree flowering inten- 
sity (Nishikawa et al, 2007; Shalom et al., 2012). In buds of 
ON-Crop trees, CIFT2 mRNA levels were relatively low, and 
remained unchanged. In contrast, CIFT2 mRNA was ~ 15- 
fold higher in buds of OFF-Crop trees. One week following 
de-fruiting, the expression of the gene in the buds of the DEF 
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Fig. 1. Fruit load affects flowering intensity the following year. Vegetative 
shoots, generative inflorescences containing only flower buds, and mixed 
inflorescences containing flower buds and leaves in various ratios were 
counted during flowering peak in trees which carried a heavy yield (ON) 
or light yield (OFF), and in de-fruited trees (DEF). The numbers are mean 
values of three independent biological replicates ±SE. Different letters 
represent a significant difference (P<0.05). 
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Fig. 2. Fruit removal alters the expression of flowering control genes in 
buds. The mRNA levels (RU, relative units) of the indicated genes were 
determined in ON-Crop (ON), OFF-Crop (OFF), and de-fruited (DEF) trees 
at the indicated weeks following de-fruiting. The numbers are mean values 
of three independent biological replicates ±SE. Different letters represent a 
significant difference (P<0.05) between the states at the same time point. 

trees was similar to that of OFF-Crop buds, and it remained 
at this level during the entire test period. The mRNA levels of 
LFY were similar in buds of ON- and OFF-Crop trees, but 
it increased 3-fold within 1 week in DEF buds and returned 
to its basal level after 4 weeks. It was previously shown that 
/w7?/56-regulated SPL5 displayed elevated mRNA levels in 
OFF-Crop buds; thus, it may act as a positive inducer of 
flowering in Citrus trees (Shalom et al, 2012). As expected, 
SPL5 mRNA levels were ~14-fold higher in buds of OFF- 
Crop trees as compared with ON-Crop buds. DEF buds 
displayed increased mRNA as compared with OFF-Crop 
buds within 2 weeks of de-fruiting. The expression of genes 
associated with trehalose and flavonoid metabolism was ele- 
vated in ON- and OFF-Crop buds (Shalom et al., 2012). As 
expected, the mRNA levels of trehalose phosphate synthase 



and trehalose phosphate phosphatase were reduced by 2- and 
4-fold, respectively following de-fruiting (Supplementary Fig. 
SI at JXB online). The mRNA levels of flavonoid biosynthe- 
sis genes, UF3GT, 4CL, CHS, and CHI, were induced by the 
treatment. De-fruiting resulted in a reduction of C4H mRNA 
levels in the ON-Crop trees. 

Fruit removal induces rapid changes in the bud global 
gene expression 

The results indicated that the transition of the bud from an 
ON to an OFF state took place relatively quickly. In order 
to analyse the metabolic and regulatory pathways playing 
a role in this transition, changes in global gene expression 
were analysed in buds before de-fruiting (Time 0) and 1 
(Time 1), 2 (Time 2), and 4 (Time 4) weeks after de-fruiting 
(Fig. 3A). In addition, buds of OFF-Crop trees were also 
analysed at Time 0 and Time 4. The minimum number of 
reads per sample was ~15 million and the maximum was ~40 
million, indicating a deep and satisfactory coverage of the 
existing transcripts (Supplementary Table S3 at JXB online). 
Overall, the number of transcripts in all libraries was ~14 
400. Hierarchical cluster analysis (Fig. 3B, left) and 2D PCA 
(Fig. 3B, right) were performed based on the log-FPKM 
values of each sample. Results of both analyses showed 
that during all time points, including Time 1 , the transcript 
profiles of DEF buds were more closely related to those of 
OFF-Crop buds than to those of ON-Crop buds, thus sup- 
porting the notion that the transition from an ON bud to an 
OFF bud was relatively quick following de-fruiting. In order 
to analyse the metabolic and regulatory pathways mediat- 
ing the ON bud to OFF bud transition, two major compari- 
sons were made (the ratios for all possible comparisons in 
the experiment are presented in Supplementary File 2 at JXB 
online). First, Time 4 included the three fruit load states, ON, 
OFF, and DEF. Therefore, the genes that were up- or down- 
regulated (P<0.05) in OFF or DEF buds in comparison with 
ON-Crop buds were identified (Fig. 4A; Supplementary File 
2). Overall, 997 genes in OFF-Crop buds and 797 genes in 
DEF buds were down-regulated relative to ON-Crop buds at 
Time 4, with 615 genes common to the two groups (OFF and 
DEF). Overall, 959 genes in OFF-Crop buds and 920 genes 
in DEF buds were up-regulated relative to ON-Crop buds 
at Time 4, with 564 genes common to the two groups (OFF 
and DEF) (Fig. 4A; Supplementary File 2). The second com- 
parison aimed at identifying genes which showed alternation 
in their expression during the course of the experiment and 
also common or different pathways altered developmentally. 
For that, genes were clustered according to their expression 
patterns relative to Time 0. For ON-Crop and DEF buds, 
four clusters were identified (Fig. 4B; Supplementary File 
2). For OFF-Crop buds, only two time points were ana- 
lysed, and altered genes were either reduced or repressed 
(Supplementary File 2) 

In order to identify common and unique metabolic and 
regulatory pathways in DEF and OFF-Crop buds which 
were altered in comparison with ON-Crop buds at Time 4, 
up-regulated genes (564+395 and 564+356, Fig. 4A) and 



Fruit load and citrus flowering | 3033 



Tree status\Samples collection Time 0 Time 1 Time 2 Time 4 

ON-Crop(ON) X X X X 

De-fruited (DEF) XXX 

Off-Crop (Off) X X 




Fig. 3. Global gene expression in buds of de-fruited (DEF) trees is similar to that of buds of OFF-Crop (OFF) trees. The experimental design and sample 
collection (A). Hierarchical clustering heatmap (left panel) and 2D principle component analysis (right panel) plots were generated using the \og-FPKM 
value of each gene for all samples, as indicated in A (B). (This figure is available in colour at JXB online.) 
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Fig. 4. Venn diagrams of differentially expressed genes and clustering analysis of developmentally altered genes. The number of genes in buds of OFF-Crop 
(OFF) and de-fruited (DEF) trees which were significantly (P<0.05) up-regulated (left) or down-regulated (right) compared with ON-Crop buds at Time 4 (A). 
Clustering analysis of genes which were significantly (P<0.05) altered developmentally during the course of the experiment in comparison with Time 0 as 
generated by Expender (http://acgt.cs.tau.ac.il/expander/overview.html) using the Click Algorithm (B). All genes are listed in Supplementary File 2 at JXB online. 
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down-regulated genes (615+382 and 615+182, Fig. 4A) were 
GO annotated and further analysed using SEA (AgriGO). 
Table 1 presents these common processes. However, in a few 
cases, unique GO terms were included under the same gen- 
eral process either because of redundancy or because they 
were part of one process. The analysis indicated that the 
great majority of common processes which were enriched in 
DEF/OFF-Crop buds were related to light sensing (such as 
response to light stimulus/intensity and response to far red/ 
blue and red light), photosynthesis (such as dark and light 
reactions and photosynthetic electron transport chain), chlo- 
roplast reorganization (such as plastid localization, chloro- 
plast relocation and organization), response to carbohydrate 
stimulus (i.e. sucrose and disaccharide), ion homeostasis 
(such as cation homeostasis and transport, proton trans- 
port), and the pentose-phosphate cycle. Fewer genes were 
down-regulated in DEF/OFF-Crop buds and belonged to 
several secondary metabolic pathways such as the terpe- 
noid and phenylpropanoid metabolic process, and oxidative 
reduction. Major unique processes induced in OFF-Crop 
buds included starch biosynthesis, carbohydrate metabolism, 
glycoside and glycosinolate metabolism, and water homeo- 
stasis (Supplementary Table S4 at JXB online). Considerably 
more unique processes were induced in DEF buds as com- 
pared with OFF-Crop buds (Supplementary Table S4). 
They included regulation of peptidase activity, regulation of 
dephosphorylation, and salicylic acid metabolism. Down- 
regulated unique processes in DEF buds included amino 
acid/amine metabolism and aromatic compound metabo- 
lism (Supplementary Table S4). Responses to biotic stress 
were either induced or repressed in DEF buds. Unique genes 
down-regulated in OFF-Crop buds could not be GO anno- 
tated, due to their low number. 

GO annotation and SEA were also performed to develop- 
mentally altered genes compared with Time 0 in ON-Crop, 
OFF-Crop, and DEF buds (Fig. 4B). Significantly altered bio- 
logical processes (FDR <0.05) could be identified for genes of 
clusters 3 and 4 (DEF buds), cluster 2 (ON-Crop buds), and 
up/down regulated genes of OFF-Crop buds. No significant 
biological processes could be identified even under FDR <0.1 
in genes of the other clusters. Cluster 4 of DEF buds included 
processes involved in light responses and photosynthesis 
(Supplementary Table S5 at JXB online). A considerable induc- 
tion was also detected in genes of respiratory burst and in those 
involved in responses to abiotic stresses. Genes up-regulated in 
OFF-Crop buds from Time 0 to Time 4 were annotated into 
light responses, photosynthesis, and response to carbohydrate 
stimulus (Supplementary Table S6). Down-regulated genes in 
OFF-Crop buds from Time 0 to Time 4 included responses 
to jasmonic acid (JA) and ethylene (Supplementary Table S7). 
Other processes included amino acid and amine metabolism, 
and responses to abiotic stresses. Cluster 3 of DEF buds included 
104 genes that could be annotated into different biological pro- 
cesses, a few of them related to development and morphogen- 
esis, such as trichoblast and root differentiation, floral organ 
development, and shoot morphogenesis (Supplementary Table 
S8). Cluster 2 of ON-Crop buds was enriched in genes associ- 
ated with lipid transport (Supplementary Table S9). 



An additional analysis for developmentally regulated genes 
was aimed at identifying common genes in clusters showing a 
similar expression pattern in OFF-Crop and DEF buds and 
opposite patterns in ON-Crop buds. The first comparison 
included cluster 4 of DEF buds, up-regulated genes of OFF- 
Crop buds, and cluster 4 of ON-Crop buds (Supplementary 
Fig. S2A at JXB online). However, only one unknown gene 
was common among the three states. The second com- 
parison included cluster 1 of DEF buds, down-regulated 
genes of OFF-Crop buds, and cluster 3 of ON-Crop buds 
(Supplementary Fig. S2B). Only the expression of three 
genes was common; two of them were homologous to terpene 
(nerolidol) syntheses (Cs2g07240 and Cs2g07250). 

Photosynthetic genes and proteins are up-regulated in 
response to fruit removal 

Consistent with a previous report (Shalom et al, 2012), the 
above-described results and MapMan analysis of differentially 
expressed genes in DEF versus ON-Crop buds at Time 4 showed 
that genes associated with light reactions, the Calvin-Benson 
cycle, and to a lesser extent photorespiration were induced 
(Supplementary Fig. S3 at JXB online). Fold changes of 41 
photosynthetic genes in DEF and OFF-Crop buds relative to 
ON-Crop buds are shown in Fig. 5A. Out of 41 differentially 
expressed genes, 38 were up-regulated while only three were 
down-regulated at at least one time point throughout the experi- 
ment. Validation of the above results by qPCR analyses was per- 
formed for seven genes, Ribulose bisphosphate carboxylase-small 
subunit (RbcS), Light-harvesting chlorophyll B-binding protein 
(LHCB3), Phosphoribulokinase (PRK), Photosystem II reaction 
centre PSB28 (PSB28), Photosystem I subunit D (PSAD), Serine 
hydroxymethyltransferase (SHM), and Ferredoxin (Fd) (Fig. 5B). 
The mRNA levels of all these genes were significantly higher in 
buds of OFF trees as compared with ON trees at Time 4. As 
expected, gene expression increased in DEF buds by 2.6- to 
5.2-fold within 1 week of fruit removal. Furthermore, western 
blot analyses using specific antibodies raised against RbcS, Fd, 
LHCII, and PsbA showed that their protein level was higher 
in OFF-Crop buds relative to ON-Crop buds, and increased 
in DEF buds (at Time 4) following fruit removal (Fig. 5C). 
However, RbcL protein levels remained unchanged in DEF and 
OFF-Crop buds relative to ON-Crop buds. 

Fruit removal-induced changes in ABA-metabolizing 
genes 

Results of genomic analysis showed that genes homologous 
to 9-cis-epoxycarotenoid dioxygenase (NCED), coding for this 
rate-limiting enzyme of ABA biosynthesis (Supplementary 
Fig. S5 at JXB online), were higher in buds of OFF-Crop 
buds relative to ON-Crop buds (Fig. 6A). In Arabidopsis, the 
NCED gene family comprises nine members, and the roles 
of NCED2, 3, 5, and 6 in ABA biosynthesis were demon- 
strated (Tan et al, 2003). The Citrus genome contained nine 
highly homologous genes, and changes in three of them 
were detected in genomic analysis: Cs5gI4370.1, homolo- 
gous to NCED3, Cs5gl4370.1, homologous to NCED1, and 
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Cs7gl4820.1, homologous to NCED4. Among these three 
genes, NCED3-like induction was increased by ~4-fold in 
OFF-Crop buds relative to ON-Crop buds at Time 0, and 
was ~2-fold higher at Time 4. Fruit removal induced a 3-fold 
increase at Time 1 and its mRNA level remained higher rela- 
tive to ON-Crop buds throughout the experiment. Induction 
of NCEDl-like was seen only in OFF-Crop buds at Time 4, 
while NCED4-like was induced in DEF buds at Time 1 and 
4 and in OFF-Crop buds at Time 4. Among the three genes, 
NCED3 is considered the major enzyme catalysing the rate- 
limiting step in ABA biosynthesis, and the mRNA levels of 
its Citrus counterpart were higher in OFF and DEF buds 
throughout the experiment. qPCR validation showed that 
Cs5gl4370.1 mRNA levels in OFF-Crop buds were signifi- 
cantly higher relative to ON-Crop buds, by ~4.5-fold at Time 
0 and by 2.1 -fold at Time 4 (Fig. 6B). The mRNA levels in 
DEF buds increased relative to ON-Crop buds by 3.2-fold 
at Time 1 and remained higher throughout the experiment. 
The transcriptomic data showed that the expression of a gene 
homologous to PYR1, a component of the ABA receptor, 
was reduced (Supplementary Fig. S6). 

The levels of ABA and its catabolites were analysed at 
three time points following de-fruiting (Fig. 7). In addition 
to ABA and its isomer, trans-ABA (t-ABA), four catabolites 
were detected in the buds, 7'OH-ABA, ABA glucose ester 
(ABAGE), phaseic acid (PA), and dihydrophaseic acid (DPA), 
with t-ABA, PA, and ABAGE showing relatively high levels. 
ABA levels in ON-Crop buds were significantly higher rela- 
tive to OFF-Crop buds, by 2.4-fold at Time 0 and by 3.3-fold 
at Time 4. In DEF buds, ABA levels decreased relative to 
ON-Crop buds by 1.3-fold at Time 1 and by 1.8-fold at Time 
4. In addition, the levels of PA and ABAGE were significantly 
higher in ON-Crop buds relative to OFF-Crop buds by 4.5- 
fold and 4.6-fold at Time 1 and by 7.7-fold and 2.9-fold at Time 
4, respectively. In DEF buds, PA and ABAGE levels decreased 
by 3.4-fold and 1.8-fold at Time 1 and by 4.1-fold and 3.5-fold 
at Time 4, respectively. The levels of DPA and t-ABA were gen- 
erally higher in ON-Crop buds relative to OFF and DEF buds, 
especially at Time 4, but the differences were not significant. 

Induction of calcium-related genes associated with 
auxin transport in DEF and OFF-Crop buds 

Among the genes which showed relatively high levels of 
expression in OFF-Crop and DEF buds relative to ON-Crop 
buds were calcium-related genes (Fig. 8A; Supplementary 
File 2 at JXB online). The expression of most of these genes 
was induced 3- to 40-fold at all time points of the experi- 
ment. Many of these genes encode proteins containing an 
EF-hand domain. PINOID (PID)-binding protein (PBPs) 
are a subgroup of EF-hand proteins. PBP1 has been shown 
to play a role in auxin polar transport in response to changes 
in calcium levels (Benjamins et a!., 2003). Four citrus PBP1- 
like genes (out of the five found in the citrus genome) were 
induced in DEF and OFF-Crop buds relative to ON-Crop 
buds at all time points tested (Fig. 8A). Phylogenic analysis of 
these genes showed remarkable homology with PBP1 and its 
closely related gene in Arabidopsis (Supplementary Fig. S8A). 
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Fig. 5. The photosynthetic machinery is up-regulated in buds of OFF-Crop (OFF) and de-fruited (DEF) trees. Fold-change (P<0.05) in the expression 
of photosynthetic genes (determined by MapMan analysis, see Supplementary Fig. S2 at JXB online) in buds of OFF and DEF trees relative to buds 
of ON-Crop (ON) trees at the indicated time points. Asterisks mark genes selected for validation by qPCR analyses, and specific genes are listed in 
Supplementary File 2 (A). Expression analysis of selected genes at the indicated weeks following de-fruiting, as determined by qPCR analyses (RU, 
relative units). The numbers are mean values of three independent biological replicates ±SE. Different letters represent a significant difference (P<0.05) 
between the states at the same time point. The lower right graph shows the linear regression between the transcriptomic and transcriptional (qPCR) 
data (B). Immunoblot analyses in two independent replicates of photosynthetic proteins (RbcL, Rubisco large subunit; RbcS, Rubisco small subunit; Fd, 
ferredoxin; LHCII, light-harvesting complex II; Psba, D1 protein) extracted from two replicates of buds of ON-Crop, (ON1 and ON2), OFF-Crop (OFF1 and 
OFF2), and de-fruited (DEF1 and DEF2) trees at Time 4 (C). The quantification of the protein signals, generated using ImageJ software, is presented in 
Supplementary Fig. S4. 



Two other induced genes, Cs9g20300.1 and Cs8g20 150.1, 
showed high homology (80% and 70%, respectively) to a cal- 
cium-dependent protein kinase (Atlg08650.1), and to a Ca 2+ / 
H + antiporter CAX3 (At3g51860), respectively. Another gene, 
Cslg21460.1, which showed homology to members of the 
NPH3 gene family from Arabidopsis also plays a role in auxin 
polar transport (Furutani et al. , 20 1 1 ; Knauer et al. , 20 1 1 ; Li 
et al, 2011; Wan et al, 2012). This family contains 33 genes in 
the Arabidopsis genome, and the Citrus genome comprises 25 
homologous genes with relatively close taxonomic relation- 
ships (Supplementary Fig. S8B). The NPH3-like gene was 
induced 25- to 40-fold in OFF-Crop and DEF buds relative to 
ON-Crop buds at all time points tested (Fig. 8A). The plausi- 
ble mechanistic relationships between calcium, PBP1, NPH3, 
and the polar subcellular localization of PIN-FORMED 
(PIN) auxin efflux carriers are schematically presented in 
Fig. 8B. The expression levels of CAX-like, NPH3-like, two 
P BP 1 -like genes, and one Ca 2+ -binding EF-hand gene were 



validated by qPCR analyses (Fig. 8C). The mRNA levels 
of all these genes were relatively low in ON-Crop buds and 
remained low throughout the experimental period. In OFF- 
Crop buds, they were significantly higher at Time 0 and Time 
4 (by factors of 5-19 and 10-28, respectively). One week after 
de-fruiting, the expression of all these genes was significantly 
increased, and within 2 weeks they attained their maximal 
levels. While the mRNA levels of CAX-like, NPH3-like, and 
Ca 2+ -binding EF-Handl-like remained high 4 weeks after de- 
fruiting, those of PBPl-likel and PBPl-like2 were reduced 
by ~2-fold and 1.5 -fold, respectively, during the same period. 

Auxin levels are significantly higher in ON-Crop buds 
and decrease following fruit removal 

Next, levels of endogenous IAA were examined in the buds 
(Fig. 9). IAA levels were significantly higher in ON-Crop 
buds relative to OFF-Crop buds by 2.9-fold at Time 0 and by 
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Fig. 6. NCED-Wke genes are induced in buds of de-fruited (DEF) trees. Fold change (P<0.05) in the expression of three NCED-like genes in buds of 
OFF-Crop (OFF) and DEF trees relative to buds of ON-Crop (ON) trees (an asterisk marks the NCED3-like gene which was validated by qPCR) (A). 
Expression of the NCED3-like gene in buds of ON, OFF, and DEF trees as determined by qPCR analyses at the indicated weeks following de-fruiting (RU, 
relative units). The numbers are mean values of three independent biological replicates ±SE (B). Different letters represent a significant difference (P<0.05) 
between the states at the same time point. 
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Fig. 7. Fruit load affects the levels ABA and its catabolites in the bud. ABA and ABA catabolites (7'OH ABA, t-ABA, frans-ABA; ABAGE, ABA glucose 
ester; PA, phaseic acid; DPA, di-hydrophaseic acid) whose metabolic relationships are schematically represented in the inner scheme were determined 
at the indicated time points in buds of ON-Crop (ON), OFF-Crop (OFF), and de-fruited (DEF) trees. Means followed by different letters are significantly 
different (P>0.05) according to Tukey-Kramer multiple comparisons tests. 
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Fig. 8. Fruit load affects the expression of Ca 2+ -related auxin polar transport and NPH3-like genes. Fold change (P<0.05) in the expression of 
Ca 2+ -related and NPH3-like genes in buds of OFF-Crop (OFF) and de-fruited (DEF) trees relative to buds of ON-Crop (ON) trees (asterisks mark 
genes selected for validation by qPCR analyses; all genes are listed in Supplementary File 2 at JXB online) (A). Expression of the indicated genes 
in buds of ON, OFFm and DEF trees as validated by qPCR at the indicated time points following de-fruiting (RU, relative units). The numbers 
represent mean values of three independent biological replicates ±SE. Different letters represent a significant difference (P<0.05) between 
the states at the same time point. (B) Linear regression between transcriptomic data and transcriptional data (qPCR analyses) is presented in 
Supplementary Fig. S7. Schematic model representing Ca 2+ and NPH3 regulation of PIN cellular localization (C). (This figure is available in colour at 
JXB online.) 
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Fig. 9. Fruit load affects auxin level in buds. Auxin was determined at 
the indicated time points in buds of ON-Crop (ON), OFF-Crop (OFF), and 
de-fruited (DEF) trees. The numbers are mean values of three independent 
biological replicates. Means followed by different letters are significantly 
different (P>0.05) according to Tukey-Kramer multiple comparisons tests. 
(This figure is available in colour at JXB online.) 



5.2-fold at Time 4. As expected, IAA levels decreased in DEF 
buds relative to ON-Crop buds, by a factor of 2.4 at Time 1, 
and remained at this level 4 weeks after de-fruiting. 

Discussion 

Intensity of return bloom is affected by de-fruiting, 
which induces relatively rapid changes in expression of 
flowering control genes and in the transcriptome 

Fruit removal has been reported to be effective in inducing 
return bloom (Monselise and Goldschmidt, 1981). However, 
since annual variation and cultivar-dependent divergence may 
affect its effectiveness (Verreynne and Lovatt, 2009; Martinez- 
Fuentes et riL, 2010; Munoz-Fambuena et al, 2011), fruit 
removal was carried out as early as the last 10 d of August. 
Indeed, the number of inflorescences and vegetative shoots, 
counted during the following spring (Fig. 1), showed that 
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those of DEF trees were similar to those of OFF-Crop 
trees, thus demonstrating the effectiveness of the treatment. 
Although it could be expected that the number of generative 
inflorescences would be higher in DEF than OFF-Crop trees 
(Verreynne and Lovatt, 2009), the actual number was identi- 
cal, probably due to the low number of fruits on OFF-Crop 
trees (~45-fold lower than in ON-Crop trees). De-fruiting 
resulted in relatively rapid changes in the expression of genes 
controlling and genes associated with trehalose and flavonoid 
biosynthesis, allowing determination of the time frame of the 
genomic analysis. 

As a rule, the expression of flowering control genes 
is induced in leaves, buds, and stems in association with 
the onset of the flowering induction period (November- 
December) in regular bearer cultivars, and in AB cultivars 
during the OFF-Crop year, while GA treatment reduced their 
expression (Nishikawa et al, 2007; Munoz-Fambuena et al, 
2011, 2012ft; Shalom et al, 2012; Goldberg-Moeller et al, 
2013). GSPL5 is an exception to this rule, probably due to its 
highly regulatory role, and it exhibited higher expression in 
buds from May until September (Shalom et al, 2012). CiFT2 
expression was usually earlier than the onset of the flowering 
induction period (Shalom et al, 2012). Taking into account 
these expression patterns and the expected year to year alter- 
nation, the mRNA levels of CiFT2 and GSPL5 were higher 
in OFF- than in ON-Crop buds. Therefore, the increase in 
the expression of CiFT2 and CiSPL5 in buds of DEF trees 
to levels similar to those in buds of OFF-Crop trees can be 
expected, and reinforces their role in return bloom. It is not 
surprising that during the time of the experiment, LFY did 
not show any difference in its expression between buds of 
ON- and OFF-Crop trees. Nevertheless, de-fruiting resulted 
in a 2-fold increase in its mRNA level, which returned to its 
basal level 4 weeks after treatment. Whether this temporary 
response has any relationship to return bloom requires fur- 
ther research. 

The differences between bud populations — those with 
a 55% chance to flower (ON-Crop) and those with a 96% 
chance (OFF-Crop and DEF buds) — did not seem to be 
very high. However, genomic analysis resulted in numerous 
differentially expressed genes (DEGs), allowing the partial 
identification of mechanisms that convert ON into OFF 
buds. Previously it was shown that the number of DEGs 
between ON- and OFF-Crop buds was considerably lower 
in September than in May (Shalom et al., 2012). However, 
the present work showed that the number of these DEGs 
was quite high. This difference could be explained by year to 
year alternations, and differences in the methodologies used. 
Based on a cut-off of 50% coverage between sequences on 
the microarray and the currently identified sequences, and on 
at least 75% identity, it is estimated that only ~30-35% of the 
current sequences are present on the microarray, supporting 
this notion. Only 40% of the auxin transport-related genes 
(Fig. 8A) were found on the microarray. Below, three identi- 
fied mechanisms, common to OFF and DEF buds, which are 
altered during the conversion of DEF buds into OFF buds, 
and might play a role in the signalling mechanism of fruit 
load are discussed. 



Induction of photosynthetic gene expression and 
protein levels in the bud following de-fruiting 

In agreement with Shalom et al. (2012), this study demon- 
strates that de-fruiting induces expression of photosynthetic 
genes in the bud. Although a recent proteomic analysis did 
not show an increase in photosynthetic proteins in OFF- 
Crop trees (Munoz-Fambuena et al, 2013), here the induced 
gene expression resulted in increased protein levels of four 
major genes. According to the C/N theory, the proteins of 
photosynthetic machinery represent the majority of leaf 
nitrogen which is directly related to photosynthetic capac- 
ity (Evans, 1989); thus, the induced levels of photosynthesis 
proteins would suggest the induction of photosynthesis in 
OFF buds, although direct evidence is missing. Although 
bud photosynthesis was never measured in fruit trees, leaf 
photosynthesis in relation to fruit load has been measured 
in previous studies. While some workers found no change in 
photosynthesis between leaves of ON- and OFF-Crop trees 
(Roper et al, 1988; Monerri et al, 2011; Nebauer et al, 
2013), others reported increased photosynthetic and C0 2 
assimilation rates in fruit-bearing as compared with non- 
fruit-bearing trees (Fujii and Kennedy, 1985; Dejong, 1986; 
Gucci et al, 1995; Palmer et al, 1997; Iglesias et al, 2002; 
Syvertsen et al, 2003; Urban et al, 2004). Vegetative growth 
is induced in buds of OFF-Crop and DEF trees (Monselise 
and Goldschmidt, 1982), suggesting that increased photosyn- 
thesis may mark the initiation of vegetative growth. That is, 
due to fruit absence, the OFF-Crop and DEF trees are heavily 
loaded with photoassimilates, suggesting that by induction of 
its photosynthetic machinery, the bud signals to stop trans- 
location of photoassimilates. The possibility that the flow of 
photoassimilates into the bud is reduced due to lower leaf 
photosynthesis in OFF-Crop trees, resulting in increased syn- 
thesis of photosynthesis proteins and higher C0 2 assimila- 
tion, cannot be excluded. 

Bud ABA is reduced in OFF-Crop trees and following 
de-fruiting compared with ON-Crop trees 

Increased expression of three NCED-like genes, in buds of 
DEF and OFF-Crop trees compared with buds of ON-Crop 
trees, suggests the induction of ABA biosynthesis. However, 
direct measurements of ABA and its catabolites showed the 
opposite trend, namely reduced levels in buds of OFF-Crop 
and DEF trees. Direct biochemical evidence demonstrated 
that NCED3 (Cs5gl4370) cleaved 9-cw-violaxanthin to form 
xanthoxin, a precursor of ABA (Kato et al. , 2006); its expres- 
sion paralleled ABA levels in the peel and during cycles of 
drought and re-watering of leaves and fruit (Rodrigo et al, 
2006; Agusti et al. , 2007). Therefore, one would expect higher 
ABA levels in OFF-Crop and DEF buds than in ON-Crop 
buds. A possible explanation of these apparently contradic- 
tory results is that the source of ABA in the ON-Crop bud 
is not within the bud itself, but external to it, and depend- 
ent on the presence of fruit. In OFF-Crop trees or follow- 
ing de-fruiting, the translocation of ABA from this source 
into the bud is blocked, at least partially, reducing the bud's 
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ABA contents, and inducing NCDE3 expression in order to 
increase endogenous ABA production. Nevertheless, the pos- 
sibility cannot be excluded that induced expression of NCED 
genes is futile, and has no physiological role. A closer look 
at ABA-responsive genes in the transcriptomic data (overall, 
eight differentially expressed genes between the states) did not 
solve this contradiction, as no common trend in their response 
was evident (data not shown). Regardless, ABA levels (and 
the expression of ABA receptor component, PYRl-like) in 
OFF-Crop buds and following de-fruiting were reduced, 
raising the question of its possible involvement in AB con- 
trol. Consistent with the present results, buds of ON-Crop 
trees have been shown to contain higher levels of ABA or its 
isomer, t-ABA, than those of OFF-Crop trees (Jones et at, 
1976; Goldschmidt, 1984). It has been suggested that elevated 
levels of ABA in ON organs may reflect a stress imposed by 
the fruit overload. Moreover, ABA might serve as an inhibi- 
tor of return bloom, since the local application of ABA to 
Citrus unshiu buds in late December inhibited bud sprouting 
and intensive flowering (Garcialuis et at, 1986). Alternatively, 
the possibility that flowering promotes ABA activity has been 
suggested, since increased ABA levels were detected in leaves 
of OFF-Crop trees and following de-fruiting of ON-Crop 
trees in association with flowering induction (Koshita et at , 
1999; Okuda, 2000). Whether ABA plays a role in AB con- 
trol, or in other processes, such as maintaining the bud in 
an inactive state (Little and Edit, 1968; Horvath et at, 2003; 
Shalom et at, 2012), requires further investigation. 

De-fruiting induces genes of calcium-dependent auxin 
polar transport 

The results showed an increase in the expression of calcium- 
related genes together with significant reduction in auxin lev- 
els in OFF-Crop buds and in buds following de-fruiting as 
compared with ON-Crop buds. Changes in the concentration 
of cytosolic free Ca 2+ ([Ca 2+ ] cyt ), mediated by ion channels, 
Ca 2+ -ATPases, and Ca 2+ /H + antiporters, form the basis of the 
Ca 2+ signalling mechanism. The CAX-type antiporters are a 
family of cytosolic low-affinity Ca 2+ /H + antiporters, which 
in Arabidopsis comprises six members. In Citrus there are 
four highly CAX homologous genes, and the expression of 
a CAX3 homologue was highly induced following de-fruit- 
ing. Transduction of Ca 2+ signals is carried out by specific 
calcium-binding proteins, containing a common structural 
motif called the 'EF-hand', a helix-loop-helix structure that 
binds a single Ca 2+ ion (Day et at, 2002). The present results 
showed a significant up-regulation of a few genes encoding 
EF-hand proteins in OFF and DEF buds compared with 
their level in ON buds. Overall, these results might suggest 
that [Ca 2+ ] cyt is affected by fruit load, although at this stage a 
plausible scenario as to the nature of the change and its cel- 
lular signature cannot be suggested. How are these changes 
related to auxin polar transport? Four of the up-regulated 
EF-hand genes show remarkable homology to the genes 
encoding PBP1 in Arabidopsis. PBP1 interacts physically 
with PID protein kinase, regulating its activity in response 
to changes in calcium levels (Benjamins et at, 2003). PID 



regulates the polarity of PIN proteins (Friml et at, 2004), 
which are known to direct auxin flow (Wisniewska et at, 
2006). NPH3-like proteins have recently been shown to affect 
PIN localization (Furutani et at, 2011; Wan et at, 2012). As 
shown here, NPH3-like genes are part of a relatively large 
gene family. Divergence of its different members occurred 
before the divergence of Arabidopsis and Citrus. The citrus 
NPH3-like gene induced in OFF and DEF buds compared 
with ON buds showed very close homology to an Arabidopsis 
gene which has not yet been subjected to detailed analysis. 

Taken together, the present results lead to the suggestion 
that higher levels of IAA in ON buds reflect their inability to 
distribute IAA efficiently via the Ca 2+ -dependent PIN-based 
polar auxin transport mechanism. In addition, efficient auxin 
removal from the bud appears to be a key component in trans- 
forming the ON bud into an OFF bud. The involvement of 
auxin in flowering inhibition following an ON-Crop year was 
recently suggested (Smith and Samach, 2013), and is based on 
the ATA hypothesis proposed by Bangerth (Bangerth, 1989, 
2006; Caaejas and Bangerth, 1997). The application of auxin 
polar transport inhibitors resulted in flowering induction in 
a number of fruit trees (Bukovac, 1968; Ben-Tal and Lavee, 
1985; Ito et at, 2001; Blaikie et at, 2004; Bangerth, 2006). 
The strong polar transport of auxin from the dominant sinks 
(i.e. the fruit or the seed), as suggested by the ATA hypothe- 
sis, preventing auxin export from the bud, would explain why 
auxin levels in OFF buds and in buds following de-fruiting 
are lower than in ON buds. 

The parallel reduction in ABA and IAA levels in the bud 
would suggest cross-talk between the ABA and IAA signal- 
ling pathways. Such cross-talk interactions were suggested in 
Arabidopsis embryo axis elongation and root development 
(Belin et at, 2009; Shkolnik-Inbar and Bar-Zvi, 2010; Wang 
et at, 201 1), but not in flowering control processes. 
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Supplementary data are available at JXB online. 
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Figure SI. De-fruiting alters the expression of trehalose- 
and flavonoid-metabolizing genes. 

Figure S2. Venn diagrams of developmentally 
regulated genes. 

Figure S3. Photosynthetic genes are induced following de- 
fruiting in the bud. 

Figure S4. Quantification of protein blot results (Fig. 5C). 

Figure S5. Schematic representation of the cleavage of 
9-cis xanthophylls to xanthoxin by 9-cw-epoxycarotenoid 
dioxygenase (NCED), a key regulated step in the biosynthesis 
of ABA in plants. 

Figure S6. Changes in the expression of the P YRl-like gene. 

Figure S7. Linear regression between transcriptomic and 
transcriptional (qPCR) data of Ca-related and NPH3-like 
genes, presented in Fig. 8A and B. 

Figure S8. Genes encoding PINOID (PID)-binding protein 
1 (PBP1) and Non-Phototropic Hypocotyl 3-like (NPH3) 
show homology in Arabidopsis and Citrus. 
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Table SI. List of primers used in this study. 

Table S2. List of genes used to design probes for nCounter 
analysis (Shalom et al., 2012). 

Table S3. Statistical data of deep sequencing analysis. 

Table S4. GO categorization of unique genes up- or down- 
regulated in buds of OFF-Crop and de-fruited (DEF) trees 
relative to buds of ON-Crop trees during Time 4. 

Table S5. GO categorization of genes of cluster 4 of 
DEF buds. 

Table S6. GO categorization of up-regulated genes of 
OFF-Crop buds. 

Table S7. GO categorization of down-regulated genes of 
OFF-Crop buds. 

Table S8. GO categorization of genes of cluster 3 of 
DEF buds. 

Table S9. GO categorization of genes of cluster 2 of 
ON-Crop buds. 
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The file includes: (i) all possible comparisons between the 
treatments; (ii) the accession numbers of genes presented 
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